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Superconducting films subject to an in-plane magnetic field exhibit a gapless superconducting 
phase. We explore the quasi-particle spectral properties of the gapless phase and comment on the 
transport properties. Of particular interest is the sensitivity of the quantum interference phenomena 
in this phase to the nature of the impurity scattering. We find that films subject to columnar 
defects exhibit a 'Berry-Robnik' symmetry which changes the fundamental properties of the system. 
Furthermore, we explore the integrity of the gapped phase. As in the magnetic impurity system, 
we show that optimal fluctuations of the random impurity potential conspire with the in-plane 
magnetic field to induce a band of localized sub-gap states. Finally, we investigate the interplay of 
the proximity effect and gapless superconductivity in thin normal metal-superconductor bi-layers. 

PACS numbers: 73.23.-b, 73.50.-h, 74., 74.76.-w, 74.25.Dw, 74.62.-c 



I. INTRODUCTION 

It is well established both theoretically and experimen- 
tally that bulk s-wave superconductivity is robust with 
respect to the addition of non-magnetic impurities (An- 
derson theorem (lj). By contrast, magnetic impurities 
affect a mechanism of pair-breaking on the system which 
suppresses and eventually destroys the superconducting 
phase. Remarkably, in such systems, the quasi-particle 
energy gap vanishes more rapidly than the order param- 
eter, admitting the existence of a gapless superconduct- 
ing phase. The (mean-field) theory of the gap suppres- 
sion was explored in a classic work by Abrikosov and 
Gor'kov Q . Since this pioneering work, it is now realized 
that this scheme applies more widely, encompassing gen- 
eral time-reversal (T) symmetry breaking perturbations 
(for a review see, e.g., Refs. |||^]). In the following, we 
will consider the general phenomenology of a disordered 
superconducting thin film subject to a homogeneous in- 
plane magnetic field. 

Generally, a bulk superconductor subject to a weak 
magnetic field of magnitude H < H C 2 enters a Meissner 
phase Q: the field is expelled from the bulk, penetrat- 
ing only a thin surface layer down to the London pene- 
tration depth. In the Meissner phase, the quasi-particle 
properties of the bulk states are not affected by the mag- 
netic field and, thus, are largely insensitive to disorder. 
However, if the dimensions of the superconductor are di- 
minished to a scale comparable to the penetration depth, 
field lines enter the sample and the superconductivity be- 
comes strongly suppressed. Here the quasi-particle prop- 
erties deviate significantly from those of the unperturbed 
system: in particular, the system can exhibit a gapless 
phase. At the level of mean-field, the properties of the 
thin film superconductor mirror those of the magnetic 
impurity system and are described by the Abrikosov- 
Gor'kov (AG) phenomenology ||. However, beyond the 
level of mean-field, phase coherence effects due to nor- 
mal disorder strongly influence the long-range properties 



of the system. In the following, we will investigate the 
properties of the gapless phase in the thin film supercon- 
ducting system focusing on two situations: in the first 
case, we will consider an arrangement of normal impuri- 
ties which are drawn at random from a 5-correlated distri- 
bution. In the second case, we will consider the supercon- 
ducting system subject to a disordered array of columnar 
defects, i.e. an impurity potential which does not depend 
on the coordinate perpendicular to the plane. The two 
arrangements are depicted schematically in Fig. [l]. 

When applied to a superconductor, a magnetic field 
has two effects: firstly, it induces a Zeeman splitting 
and, secondly, it couples to the orbital motion of the 
electrons. Both have the effect of suppressing supercon- 
ductivity. Whereas, in a bulk system, the orbital effect 
usually dominates, in very thin films the opposite situa- 
tion arises. The crossover can be estimated in the follow- 
ing way ||: the critical magnetic field associated with 
the orbital effect is roughly determined by the condition 
that the flux threading an area spanned by the coherence 
length is of the order of one flux quantum, H^ 2 ~ <fio, 
where £ = \J D I (2A) is the coherence length in the dirty 
system with order parameter A and diffusion constant D. 
Now, if the width of the film d becomes smaller than the 
coherence length, d <C £, this has to be replaced by the 



where H^ 2 represents the in-plane 



condition H^d 
field. I.e., the orbital critical field increases. The critical 
field Hz associated with the Zeeman splitting is inde- 
pendent of the width of the system. Hz is obtained from 
the condition that the energy splitting between up(|)- 
and down(J,)-spins is roughly of the size of the order pa- 
rameter, gL/J'BHz — A, where gh is the Lande g-factor 
and (is — e/(2m) the Bohr magneton. Comparing the 
two, one concludes that the orbital effect is dominant in 
suppressing superconductivity when 

d> d c = g L — , 
where Af is the Fermi wavelength and t the elastic mean 
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free path. 

In the following, we restrict attention to systems where 
d > d c and the Zeeman splitting can be neglected. To 
be concrete, we will consider a thin film system where 
the relevant length scales are arranged in the following 
hierarchy: 



a f < [d, i] < e 



(1) 



The inequality Ap <C I defines the quasi-classical regime 
while £ <gc ^ is the condition for the dirty limit. Finally, 
Af -C d implies that the subband splitting due to size 
quantization is small and many subbands are occupied. 
Thus, from the point of view of the normal electrons, the 
system is effectively three-dimensional. 




FIG. 1. Schematic picture of a thin film with 5-correlated 
disorder (left) and columnar defects (right). 



Before turning to the theoretical analysis, let us first 
briefly summarize the main results of this investigation. 
At the level of mean-field, the superconducting system 
with an in-plane magnetic field is described by the AG 
theory 0], independent of the nature of the normal dis- 
order potential (diffusive or columnar). The parameter 
governing the suppression of the quasi-particle energy 
gap is set by 
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where is the flux through an area perpendicular to 
the field spanned by the width of the film d and the su- 
perconducting coherence length £. As found for the mag- 
netic impurity system within the gapped phase, the 
sharp edge in the quasi-particle density of states (DoS) 
predicted by mean-field is softened by the nucleation of 
localized tail states which extend into the sub-gap re- 
gion. In this sense, the mean-field gap edge becomes 
a mobility edge separating a region of bulk (weakly lo- 
calized) quasi-particle states from strongly localized tail 
states. In the following, we will demonstrate that these 
sub-gap states are induced by optimal fluctuations of the 
random impurity potential and are reflected as inhomo- 
geneous instanton field configurations of the mean-field 
equations. In the vicinity of the gap edge -E ga p, the en- 
ergy scaling of these tails is universal |7| , depending only 
on the distance from E gap , the dimcnsionless parameter 
£ and dimensionality In 2d, one obtains the result, 



< -Egap) ~ exp 



(2) 



non-perturbative in the inverse dimensionless conduc- 
tance l/(v a D) : where v a is the DoS of the system when 
in the normal state. Here a 2 is a known dimensionless 
function of the control parameter £. 

At the mean-field level, the choice of disorder does not 
affect the results qualitatively. However, within the gap- 
less phase, the low-energy physics depends sensitively on 
the character of the impurity potential and, in particular, 
on whether or not it respects inversion symmetry. Within 
the diffusive regime I/r > e > D/L 2 , the quasi-particle 
states exhibit a localized or "spin-insulator" phase with 
a localization length £i oc which depends on the funda- 
mental symmetry of the system [p|-[lO|. However, it is 
within the ergodic regime e < max\Dj£,f oc , D j L 2 ] , where 
the influence of the inversion symmetry is most clearly 
exposed. Here the quasi-particle properties become uni- 
versal, dependent only on the fundamental symmetry 
class |lf],[l2|]. A complete classification scheme of ten 
symmetry classes 0, corresponding to the ten large 
families of symmetric spaces identified by Cartan, ex- 
ists; four of these symmetry classes are relevant in the 
context of superconductity jl2| . Possible symmetries are 
time-reversal, spin-rotation, particle-hole and chiral sym- 
metry. Now s-wave superconductors possess particle-hole 
symmetry but, generically, do not exhibit chiral symme- 
try. Furthermore, they may or may not possess time- 
reversal and/or spin-rotation symmetry. Here, as one 
might expect, one finds that the diffusive film in a mag- 
netic field belongs to 'symmetry class C (according to 
Cartan's notation) which corresponds to spin-rotation 
symmetry, but broken time-reversal symmetry. By con- 
trast - despite the presence of the magnetic field - the 
film with columnar defects manifests a hidden symmetry 
placing the superconductor in the higher symmetry class 
CI which is usually characteristic of systems possessing 
time-reversal invariance. The distinct symmetry classes 
manifest themselves in the energy dependence of the den- 
sity of states for e — > 0. Thus, although superficially 
the suppression of superconductivity does not respond 
to geometrical symmetries, manifestations of the Bcrry- 
Robnik symmetry effect can be observed in the gapless 
phase. The distinct low-energy behavior predicted above 
is confirmed by numerics. 

Finally, to complete our discussion, we turn to the con- 
sideration of the influence of an in-plane magnetic field on 
the spectral properties of a normal metal-superconductor 
(NS) bi-layer. Here the interplay of macroscopic quan- 
tum coherence phenomena in the superconductor and the 
mesoscopic metal leads to interesting effects |Q . Gener- 
ally, a superconductor brought into contact with a nor- 
mal system tends to impact aspects of its superconduct- 
ing character onto the latter. At the origin of this phe- 
nomenon is the Andreev reflection Jl4j of electrons at 
the interface: an electron with energy smaller than A 
may be retro-reflected off the boundary as a hole, and 
a Cooper pair added to the condensate in the supercon- 
ductor. The coherent superposition of the incident elec- 
tron with the reflected hole leads to a non- vanishing pair 
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amplitude within the normal region. In particular, the 
density of states within the normal region may develop a 
gap - whose size depends on the coupling 7 between the 
two systems. Here we consider the case of weak coupling, 
where the induced gap Eg^p = 7 (< A). A parallel mag- 
netic field counteracts this phenomenon by suppressing 
the energy gap even within the superconductor. Below, 
we will show that the field necessary to suppress the in- 
duced gap in the normal region is much smaller than 
the field required to drive the superconductor into the 
gapless phase. However, the coupling to the normal re- 
gion also induces a non-vanishing DoS at low energies 
E g aJ < e < A in the superconductor. Here we consider 
two different geometries, namely a planar NS bi-layer as 
well as an NS cylinder with the magnetic field directed 
along the axis. In the second system the enclosed flux 
leads to a periodic modulation of the field effect as first 
observed in an experiment by Little and Parks |D| . Only 
for small enough cylinders the system reaches the gapless 
phase in certain ranges of the magnetic field. By further 
decreasing the radius, superconductivity is completely 
suppressed around half-integer flux quanta threading the 
cylinder, as predicted by de Gennes fl6|l a nd very recently 
verified experimentally by Liu et al. fllTf! . 

The paper is organized as follows: in Sec. |n|, we first 
study the suppression of superconductivit y du e to the 
parallel field at the mean-field level. In Sec. |II B| we show 
how the results for the sub-gap tail states, recently estab- 
lished for the magnetic impurity system, translate to thin 



films in parallel fields. Subsequently, in Sees. [II and IV 



we search for manifestations of the Berry-Robnik symme- 
try effect on the properties of the superconducting phase. 
In Sec. [v|, we discuss the properties of thin disordered NS 
bi-layers. Finally, we conclude in Sec. [Vl|. 



II. DIFFUSIVE FILM 

Before turning to the columnar defect system, let us 
begin with a discussion of the most generic case: a thin 
film with just a 'normal' ^-correlated white noise disorder 
potential. Here we are interested in the limit, where - in 
addition to the conditions (|l|) specified above - £ <C d 
which implies diffusive motion in all three directions. 

Building on the field theory approach to the study 
of weakly disordered systems JT^-pO| (for a review see, 
e.g., Ref. 21|), the extension to the consideration of 
disordered superconducting systems follows straightfor- 
wardly ]22}-|2~4 ] . Therefore, here we will only briefly sum- 
marize the main elements in the construction of the field 
theory of the disordered superconductor in the framework 
of the non- linear cr-model (NLcrM). Using this formula- 
tion, we will thereafter investigate the response of the 
superconducting film to an in-plane magnetic field. 

The superconducting system is described by the 
Gor'kov Hamiltonian 

'Ho A 



with Hq = p 2 / (2m) — ep + V(r). Here, the matrix struc- 
ture refers to the particle/hole (ph) space, £f denotes the 
Fermi energy, and V(r) represents the quenched random 
impurity potential whose distribution is characterized by 
the mean scattering time r. We consider the potential 
to be drawn from a Gaussian white noise distribution, 
(V{r)V(r')) v = (27r^ r)- 1 ( 5( 3 )(r - r'), where v is the 
average DoS of the normal system. 

The starting point of the analysis is the functional field 
integral for the generating functional 



where \& represent 8-component supervector fields, in- 
corporating a boson/fermion space as well as the par- 
ticle/hole and a charge conjugation space: to prop- 
erly account for all channels of quantum interference, 
it is standard practice to affect a doubling of the field 
space to accommodate the particle-hole symmetry H — 
—<T2 n H T cr™i where tr™ is a Pauli matrix in particle/hole 
space. This additional space is referred to as the charge 
conjugation sector. The introduction of commuting 
(bosonic) as well as anti-commuting (fermionic) fields en- 
sures the normalization of the generating functional in 
the absence of sources, Z[Q] = 1. Following Ref. |J, the 
superfields are not independent, but obey the condition 
* = (icr^ H (g)?7*) T , where 77 = of c ® E BB - ia% c g! E FF 
(and E BB = diag(l,0) BF , E FF = diag (0, 1) 

bf are projec- 
tors onto the boson-boson and fermion-fermion block, re- 
spectively). Taking our notation from Ref. the action 
assumes the canonical form 



dr ¥(e - a§° 



H 



A* 



-H r , 



where a^ c is a Pauli matrix in charge conjugation space, 
e~ = e — iO, and H represents the Gor'kov Hamiltonian 
introduced above. 

To explore the low-energy properties of the supercon- 
ducting system, after ensemble averaging over the ran- 
dom impurity distribution, the functional integral over 
the supervector fields 'J can be traded for an integral 
involving supermatrix fields Q. Physically, the fields Q, 
which vary slowly on the scale of the mean-free path £, 
describe the soft modes of density relaxation — the diffu- 
sion modes. In the quasi-classical limit (ep 3> 1/t), the 
action for Q is dominated by the saddle-point field con- 
figuration. In the dirty limit (1/r ^> A), the saddle-point 
equation admits the solution Qo — er™®^ . However, in 
the limit e — > 0, the saddle-point is not unique but spans 
an entire manifold parameterized by the unitary transfor- 
mations Q = TQoT^ 1 . Taking into account slow spatial 
and temporal fluctuations of the fields, the low-energy 
long-range properties of the weakly disordered supercon- 
ducting system are described by an action of non-linear 
CT-model type pfS§ 
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S[Q] 



TTVq 
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J dr Str [D{dQf - 4(ie-er™ ® cr£ 



(3) 



where the supermatrix field obeys the non-linear con- 
straint Q 2 = 1. Furthermore, the matrices Q obey the 
the symmetry relation Q = a™ ® riQ T (cr™ ® 7 7) T j re ~ 
fleeting the symmetry properties of the dyadic product 
\& ® ^cr™. Finally, due to gauge invariance properties of 
the action, the incorporation of a magnetic field amounts 
to replacing the derivatives in Eq. (|3|) by their covari- 
ant counterpart d = d — iA[a^ H , . ]. This completes the 
formulation of the disordered superconducting system as 
a functional field integral involving the supersymmetric 
NLoM action. Our interest here is in the thermodynamic 
DoS obtained as 



7T- 1 / £$>[G(r,r;0] 



the fact that, in a superconductor, the vector potential 
is associated with a supercurrent j s = n s A/m, where n s 
is the density of Cooper pairs. The first condition tells 
us that no net current is generated while the second con- 
dition does not allow a supercurrent to flow through the 
superconductor-vacuum boundary. Thus, when integrat- 
ing along z, we have fixed the gauge, i.e. the resulting ac- 
tion is not gauge invariant. Therefore, the magnetic field 
does not appear within a covariant derivative, but as an 
additional diamagnetic term ~ K[a^ H ,Q] 2 . This distin- 
guishes the "thick" film, d^> Xp, from the single-channel 
case (i.e., the truly two-dimensional system), where the 
magnetic field can be gauged out and, thus, has no influ- 
ence. 



Making use of the generating functional, it is straight- 
forward to show that the impurity averaged DoS can be 
obtained from the identity 



(K e )> 
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dr 

V 



n(Str [a 



Q ' 



(4) 



-S[Q] 



where (■ ■ -)q = J DQ ■■■ e~ 

Thus, the starting point of our analysis is the conven- 
tional NLctM, Eq. (^), for a three-dimensional supercon- 
ducting system subject to a magnetic field (2^^], where 
the vector potential reads A = —Hze y . The typical scale 
of variation of the Q-fields is set by the coherence length 
£. Therefore, since d -C £, the matrices Q are constant 
along the z-direction, allowing integration along z to be 
performed explicitly. Making use of the identities 

d/2 d/2 

I / *A = 0, \ J dzA 2 = ±(Hd) 2 , 



-d/2 

one finds that 



-d/2 



s = - 



d 2 r Str 



D(dQ) 2 --[a : 



Qf 



(5) 



where K = D(Hd) 2 /6 and vu is the DoS of the two- 
dimensional system, i.e. here u a = v^d. In this planar 
geometry, the paramagnetic term, which is crucial for 
the Meissner effect, vanishes from the action. 

Note that the choice of gauge is important here: the 
physical gauge to choose is the London gauge, V • A = 
and A z (±d/2) = 0. Both conditions are fulfilled by 
A = —Hze y . The above requirements originate form 



A. Mean-field analysis 

As usual, at the mean-field level, the density of states 
is obtained by subjecting the action (|J) to a saddle point 
analysis. Varying the action with respect to Q, subject 
to the non- linear constraint Q 2 = 1, one obtains the 
saddle-point equation 



Dd{QdQ)--[a™Qa™,Q] 



,Q] = o. 



(6) 



The latter can be identified as the Usadel equation ]25| ] 
for the average quasi-classical Green function, supple- 
mented by an additional term due to the parallel mag- 
netic field. With the Ansatz Q = cosh^tr™ ® af ] + 



i sinh 9 cr™ , where 9 = 



j-Eff, one obtains 



Dd 2 9 - 2i(e sinh 9 - A cosh 9) - k sinh(2(9) =0. (7) 
Taking 9 to be homogeneous in space, and defining 



e KCOsh#, 

2 



A 



A + -Ksinh( 



the equations for e and A take the BCS form and, thus, 
admit the diagonal solution e/ A = coth 9. Then, in terms 
of the 'bare' parameters e and A, the saddle-point equa- 
tion (Q) can be brought to the conventional AG form 



e 

A 



u 1 - C 



(8) 



where u = coth0. Following for example Ref. g, by 
cxtrcmizing this equation and defining £ = k/A, one ob- 
tains the quasi-particle energy gap 
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E„ 



A(l-C 



2/3\3/2 



When combined with the self-consistency equation for 
the order parameter [|1 



A = 



ngvu 



^ ./n 



where we have shifted to Matsubara frequencies, ie~ — > 
e n = (2n + 1)tt/(3, one finds that, at T = 0, super- 
conductivity is destroyed when k\ — Ao/2, where Ao 
is the order parameter in the absence of a magnetic 
field. The gapless phase sets in at the smaller value 
Kg = Ao exp[— 7r/4]. In the following, we will be able 
to determine the order parameter self-consistently, and 
A will be understood as the renormalized order parame- 
ter, even if not stated explicitly. The suppression of the 
energy gap -E gap and the order parameter A are shown 
in Fig. 0. 



E /A n 

gap 




FIG. 2. Dependence of the energy gap i? gap and the order 
parameter A on k. Gapless sueprconductivity occurs in the 
narrow range Ao exp[— 7r/4] < k < Ao/2. 



Finally, in the mean-field approximation, the DoS 
takes the form 



1/(6) 



v u f d 



V 

v u K[cosh0] 



KStr[of ®t7f ®uf Q(t)] (9) 



Inserting the solution determined by Eq. (|S|), this yields 
the characteristic AG DoS, exhibiting a square-root edge 
in the gapped phase. 



by magnetic impurities have shown that the hard edge 
predicted by the AG theory is untenable [^6|,0. In par- 
ticular, in the diffusive system, it has been shown that 
optimal fluctuations of the random impurity potential 
lead to the nucleation of localized sub-gap quasi-particle 
states which soften the gap edge (?]]. 

Within the framework of the supersymmetric field the- 
ory, the precise mechanism of sub-gap state formation 
has been elucidated in Ref. jj. There, it was shown 
that, in addition to the homogeneous AG solution, the 
saddle-point equation (]?]) accommodates a degenerate 
manifold of spatially inhomogeneous instanton or bounce 
solutions. Referring to Ref. [j7j for details, to exponential 
accuracy, the sub-gap DoS takes the form 



v(e) ~ exp 



-c 2 VnD(' 



-2/3, 



i-C 



2/3\ 



-1/2 -^gap' 



(10) 



where C2 is a numerical constant. Within the same the- 
ory, one finds the sub-gap states are confined to droplets 
of size 



A 



ap 



1/12 



E B 



A 



-1/4 



diverging upon approaching the gap edge. The mean- 
field gap edge -Egap therefore assumes the significance of 
a mobility edge separating localized sub-gap states from 
extended bulk states. 

Being confined to a region of size rd rop ^ £ ^> t ^> Xf , 
it is evident that the physical mechanism of sub-gap 
state formation is quasi-classical in origin, relying on 
optimal fluctuations of the random impurity potential. 
Moreover, in contrast to Lifshitz semiconductor band- 
tail states |27| , each droplet leads to the nucleation of an 
entire band of localized states. Qualitatively, the phys- 
ical mechanism of sub-gap state formation is connected 
to mesoscopic fluctuations in the phase sensitivity of the 
electron wavefunction. In regions where the sensitivity 
is high, the impact of a time-reversal symmetry breaking 
such as that imposed by the magnetic impurities or paral- 
lel magnetic field is stronger, and the effective scattering 
rate C is enhanced. 

Note that, instead of applying a magnetic field, one 
might consider driving a supercurrent through the sys- 
tem. As emphasized in Ref. |3j , this leads to the same AG 
mean-field results. However, in this case, due to the pres- 
ence of a paramagnetic term in the action, the mean-field 
gap edge is robust, i.e. no sub-gap states are generated. 



B. Sub-gap states 

According to the mean-field AG theory, the supercon- 
ducting system is expected to exhibit a gapped phase 
over a wide range of parameters with the DoS varying 
as a square-root of the energy difference e — E gap in the 
vicinity of the gap edge E gap . Recently, studies in the 
related context of a superconducting system perturbed 



III. COLUMNAR DEFECTS 

In the normal system, previous studies [p8[ have shown 
that the properties of a thin film in an in-plane magnetic 
field depend sensitively on the nature of the impurity po- 
tential. Do the quasi-particle properties of the disordered 
superconductor exhibit a similar sensitivity? 
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In the normal system, z-inversion symmetry effectively 
compensates for the time-reversal symmetry breaking of 
the magnetic field, driving the system into the orthog- 
onal symmetry class p9| . As such, one might expect 
that the inversion symmetry restores the validity of the 
Anderson theorem, rendering the superconducting phase 
insensitive to the disorder potential. In fact, constraints 
imposed by self-consistency of the order parameter pre- 
vents the existence of a mechanism that could cancel the 
effect of the magnetic field. The reason is that here we are 
dealing with an interacting problem: with the formation 
of Cooper pairs, it is not possible to replace time- reversal 
by any other symmetry. Indeed, at the level of the mean- 
field AG theory, the influence of columnar defects is in- 
distinguishable from that of the diffusive scatterers - only 
the parameter £ is modified. 

However, as we will see in S ec, jivl , the Berry-Robnik 



phenomenon described in Ref. 29 is not completely in- 
effective in the superconducting phase. Taking into ac- 
count fluctuations in the gapless phase, we will show that, 
while the diffusive film belongs to the fundamental sym- 
metry class C (corresponding to a disordered supercon- 
ductor in a magnetic field), with columnar defects, the 
system belongs to the higher symmetry class CI (char- 
acteristic of the time-reversal invariant superconductor) . 
The result is a substantial modification of the low-energy 
behavior: in the diffusive film, quantum interference phe- 
nomena in the particle/hole channel induce a microgap 
structure with a DoS varying as v(e) ~ e 2 , while in the 
film with columnar defects v{e) ~ e. 

To be specific, let us consider a model of a thin film 
superconductor subject to a random (impurity) poten- 
tial which varies only along the in-plane directions. In 
the absence of a magnetic field or superconducting order 
parameter, the quasi-particle Hamiltonian can be subdi- 
vided into different subbands labeled by an index k. The 
spectral properties of each subband are described by a 



two-dimensional NLctM action of conventional type. The 
derivation of an effective low-energy action follows closely 
the normal case in Ref. pgfl . The Gor'kov Hamiltonian 
of the system now reads 



n = 



+ W(z) - V(x,y))a? + A(*)af , (11) 
V 2m / 

where d = d + iHzeya™, W is the confining potential, 
and V represents an impurity potential drawn at random 
from the white-noise 5-correlated distribution with zero 
mean, and correlation (V(r)V(r')) = (2iri>nT)~ 1 6( 2 \r — 
r'), where are 'in-plane' two-component vectors. 

Diagonalizing the z-dependent part of the problem, 
and representing Ti in the basis of the (real) cigcnfunc- 
tions {4> k }, i.e. H k k> = J dz (j) k H <fi k > where (-<9 2 /(2m) + 



W(z) 



e k 



0, the vector potential A = —Hze y as 



well as the order parameter become matrices in fc-space: 



A fefc / = -He y J dzcj) k {z)z<i>ki{z), 
Afcfc' = / dzc/)k{z)A(z)(p k >(z). 



Let us emphasize that, if the system possesses inver- 
sion symmetry z —> — z, the matrix clement A kk r differs 
from zero only if k + k' is odd; in particular, A kk = 0. 
For simplicity, here we only consider the fully symmetric 
case [p0|| . 

Under the further assumption that the subband spac- 
ing \e k — e k /\ is larger than the scattering rate jnj, one 
finds that only the diagonal components of the order pa- 
rameter are non- vanishing. Starting from the conven- 
tional superconducting 2c? NLerM action for the k sub- 
bands and turning on an in-plane magnetic field, it is 
straightforward to show that the total effective action 
assumes the form 



S = 



d 2 . 



•E 

k 



Str 



D k {dQ k f-A( 



16 (7 q 



AfcfcCT™) Qk — 2 ^ X kk > cr^ H Q k a^ H Q k 



(12) 



where X kk > = D kk , A kk , A k > k /{1 + {E kk >T) 2 ). Further- 
more, D kk i = (D k + D k >)/2 (with D k denoting the dif- 
fusion constant of subband k) and E kk i = e k — e k >. Cru- 
cially, from this result we see that there exists no linear 
coupling of Q to the vector potential - a paramagnetic 
term does not appear. 

To proceed, as before we subject the action ( |l2| ) to 
a mean- field analysis. Varying the action with respect 
to fluctuations of Q k , one obtains the modified (set of 
coupled) Usadel equations 

D k 8 (Q k dQ k ) - [ic-of <E> of - A fefc a™, Q k ] - 



I 

-J2 X ^' Wl H Qk<Jr,Qk'} =0. 

k' 

Applying the Ansatz Q k = cosh 8 k cr™ ® <jg c + i smli^ a™ 
with 8 k homogeneous in the in-plane coordinates, the 
mean-field equation assumes the form 

e~ sinh 9 k -A kk cosh 8 k - i ^ X kk > sinh (9 k + 6 k > ) = 0. (13) 

fe' 

In principle, this equation has to be solved in parallel 
with the self-consistent equation for the order parameter 

Afefc' = ^— — sin 9 k , n S kk > , (14) 



G 



where g is the effective BCS coupling constant and n is 
a Matsubara index. 

Analyzing the saddle point equation (Jl^) , it can be eas- 
ily seen that the field dependent term Jj,, Xkk' sin(0fc + 
9k') vanishes if we choose the solution 9k = ( — l) h 9 
(i.e. since Xkk 1 = for k + k' even). Thus, there seems to 
be one mode which is not affected by the magnetic field. 
However, this would imply that the order parameter, too, 
must have an alternating sign, i.e. A^ = (— l) k A. Re- 
calling the definition Akk = J (izA(z)^, this is not fea- 
sible. Thus, the solution above is ruled out |52| and, 
therefore, on the mean-field level, the symmetry mecha- 
nism is ineffective. 

A more natural choice seems to be a spatially homo- 
geneous order parameter. Unfortunately, for a general 
model, the solution of the coupled Eqs. ( |l3| ) and ( |l4| ) 
does not seem to be readily accessible analytically. How- 
ever, to gain some insight into the nature of the general 
solution, we will specialize further consideration to the 
particular case in which only the lowest two subbands 
are coupled. 

With A12 = X21 = X the equations for 9\ and 82 coin- 
cide. Therefore, setting 9 = 9\ =82, which implies that 
An = A22 = A, the mean-field equation takes the form 
reminiscent of the AG equation, 



e sinh 9 - A cosh — iX sinh (29) = 0. 



As with the diffusive film, the application of a strong in- 
plane field suppresses the order parameter and allows for 
the existence of a gapless phase. According to the AG 
theory, the superconductor enters the gapless phase when 
( = 2X/A~1. 

If -E12T <C 1, the parameter £ is of the same form 
as that found in the diffusive case, i.e. £ ~ D(Hd) 2 / A. 
In the opposite limit, £ is greatly reduced because the 
wide subband spacing restricts the motion in z-direction. 
Now, C ~ D(Hd) 2 /((Ei 2 T) 2 A), and, thus, higher mag- 
netic fields have to be applied in order to reach the gap- 
less phase. As in the diffusive case, the hard edge in the 
gapped phase is compromised due to fluctuations - see 
the discussion above - and exponentially small tails in 
the sub-gap region arise. 

More generally, for many subbands, one would ex- 
pect the same qualitative picture to hold - although Akk 
might slowly depend on k. 




J , i 

0.5 1 1.5 2 

e/|A| 

FIG. 3. Numerical results: DoS. Upon increasing the mag- 
netic field, the energy gap closes and the BCS singularity 
disappears. 



The effect of gap suppression is born out in a simple nu- 
merical simulation. Fig. || shows the quasi-particle DoS 
for a two subband tight-binding model with 2 x 20 x 20 
sites when subject to an in-plane magnetic field. The en- 
ergy is measured in units of the (unperturbed) order pa- 
rameter. The three curves correspond to different values 
of the magnetic field. Details of the result at intermediate 
fields are magnified in Fig. |^. The mean-field square-root 
edge as well as the exponentially small tails are indicated. 
Furthermore, the inset shows the linear energy depen- 
dence of the sub-gap DoS exponent, c.f. Eq. ([10|), on a 
linear-log scale: ln^(e < -E gap ) ~ -E-gap — e- 




0.1 0.2 0.3 0.4 0.5 

Energy £/ |A| 



FIG. 4. Numerical results: sub-gap DoS. The solid lines 
show the square-root edge and the exponentially small tails. 
In the inset, the same data are plotted on a linear-log scale. 
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IV. PHASE COHERENCE PROPERTIES OF THE 
GAPLESS PHASE: MASSLESS FLUCTUATIONS 
AND THE SOFT MODE ACTION 

While, at the level of mean-field, all perturbations 
(i.e. magnetic impurities as well as parallel fields in films 
with different disorder potentials) follow the same AG 
phenomenology, it is interesting to note that differences 
show up in the spectrum of soft fluctuations. The latter 
are responsible for the long-range spectral and localiza- 
tion properties of the quasi-particles in the gapless phase. 
In contrast to the magnetic impurity model j|]7j which 
belongs to symmetry class D (due to broken time-reversal 
and spin-rotation symmetry), here the soft fluctuations 
around the mean-field should be described by an effec- 
tive action belonging to symmetry class C (broken T- 



invariance) 112 
tions of Ref . 




Therefore, according to the considera- 
lO) , one expects localization of the quasi- 
particle states in the gapless phase. In fact, below we will 
see that the fluctuations are sensitive to the very nature 
of the impurity scattering. It turns out that the film with 
columnar defects is described by the higher symmetry 
class CI, implying a modified localization length |8TJlfj|l. 

To assess the low-energy properties of the system, we 
have to first identify the soft modes of the action. For fre- 
quencies e — > 0, the saddle point is not unique, but spans 
a degenerate manifold Q = TQspT" 1 with T = exp[VF] 
and {Q SP , W} = 0. The symmetries of the system impose 
certain conditions on the generators W. 



where Q s = T s o-^°T~ x and T s = exp[W s ]. 

On energy scales e < E c = D \ cosh 9\/L 2 , the system 
enters the universal zero-dimensional regime. Here the 
properties of the action are dominated by the zero spa- 
tial mode and lead to 11121 



v{e) = v(E c ) 1 



sin(27re/<5) 
27re/<5 



(16) 



where 6 = 1/ (v(E c )L 2 ). I.e. for e — > 0, the DoS vanishes 
quadratically, 



v{E c ) 3 



2 ,/e\ 2 
* '5, 



This is to be contrasted with the low-energy behavior 
of the DoS in the case of columnar defects, where the 
system possesses the ^-symmetry. 



B. Columnar defects 

Here instead of a single generator W, one has to con- 
sider a set of generators Wk ■ 

• Once again, Wk has to commute with the order 
parameter, 

w 2 \w k ] = o. 



A. Diffusive film 

The choice of generators W is dictated by the presence 
of the order parameter and the magnetic field. This leads 
to the following conditions: 

• Firstly, W has to commute with the order param- 
eter, 

{a™,W} = 0. 

• Secondly, since time-reversal symmetry is broken 
by the magnetic field, W has to fulfill the further 
restriction 

Thus, W = 1 PH ® W s . Taken together, these restric- 
tions limit field fluctuations to those belonging to sym- 
metry class C which describes superconducting systems 
with spin-rotation symmetry, but broken T-invariance. 
The corresponding integration manifold of class C is 
Osp(2|2)/Gl(l|l). Expanded in the generators, the soft 
mode action reads 

Sq s = -^p J d 2 r Str [Dcosh^SQ,) 2 + 

+4iecosh(9cr£ c Q s ], (15) 



• However, even though time-reversal symmetry is 
broken by the magnetic field, the generators do not 
have to obey [ct™, Wfe] = 0. Due to ^-symmetry, 
which causes all elements X^y with k + k' even to 
vanish, it is sufficient to require 

W k > = of 1 Wfeo| H for k + k' odd. 

I.e. one generator, take e.g. Wo, can be chosen 
'freely'. Then, the others are determined through 
the condition 

W k = (a 3 ™) fc ^oK H ) fc , 

or: W k = W if k G 2N, and W k = of H W a| H if 
k G 2N+1. 

Thus, the second condition here only imposes certain 
relations between different Wk, but does not restrict 
the structure of Wk in particle-hole space. This corre- 
sponds to the higher symmetry class CI. Now the inte- 
gration belongs to the group manifold Osp(2|2). Again 
we find a manifestation of the Berry-Robnik symmetry 
phenomenon: the low-energy properties of the gapless 
phase are determined by the symmetry class associated 
with systems possessing time-reversal invariance. 
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Taking into account these fluctuations, the correspond- 
ing soft mode action reads 



5, 



Q, 



d 2 r Str 



D k cosh 2 9 k (dQ s ) 2 



-4ie cosh.t9j.cr. 



fc°3 



(17) 



where Q s = T^fgiafTf 1 . Here T s = exp[W ], and W 
fulfills the conditions specified above. Once again, pro]: 
erties of the class CI are available in the literature 
In particular, for small energies, one obtains 



v{e) 



Tre/8 
7T f dz 



■J (z)Ji(z) 



v{E c ) 2 J z -^-'"^ ' 4 
o 

showing the DoS to vanish linearly as e — » 



e/5 + 0(e 3 ), (18) 



< 

O 
Q 




0.001 



0.010 



e/|A| 



0.100 



1.000 



FIG. 5. Numerical results: log-log plot of the low-energy 
DoS in the gapless phase for a diffusive film (open squares) 
and a film with columnar defects (open circles). 



The predicted low-energy behavior can be verified nu- 
merically. In Fig. ||, the density of states at low energies 
is compared for the two cases. On the log-log scale one 
can read off the exponent a governing the energy depen- 
dence, |e| a . At low energies, the two lines with slope 
«c = 2 and aci = 1 _ characteristic of the symmetry 
classes C and CI - fit the data for the diffusive film and 
the film with columnar defects, respectively. 

Having studied the influence of a parallel magnetic field 
on the properties of the superconducting film, we now 
turn our attention to the mean-field properties of normal 
metal-superconductor hybrid systems. 



V. NS HYBRID SYSTEMS 

The properties of thin disordered NS bi-layers have 
been studied in a recent work by Fominov and 
Feigel'man |ll|. By means of the coupled Usadel equa- 
tions for the hybrid system, they investigated the density 



of states as well as the parallel and perpendicular critical 
fields of the bi-layer as a function of the interface trans- 
parency. The asymptotics at high and low transparencies 
are accessible to analytical solutions while the results at 
intermediate transparencies were found numerically. 

Here, we consider a different aspect of the properties of 
the hybrid system, namely the interplay between gapless 
superconductivity and the proximity effect. In addition 
to the effect of the field on the individual system, in the 
NS bi-layer it also affects the coupling. As we have seen 
in the previous sections, an in-plane magnetic field grad- 
ually suppresses the gap in the single-particle DoS. On 
the other hand, in an NS structure, the proximity effect 
opens a gap in the DoS of the normal layer. Thus, we 
expect the magnetic field to weaken the proximity effect. 

For simplicity, we consider here a hybrid system each 
consisting of a single N- and S-channel (see Fig. fy. I.e., 
neglecting the finite width, the magnetic field does not 
influence the individual systems, and we can study the 
effect on the coupling alone. The coupling between the 
layers is described by a tunneling Hamiltonian H.t = 
J d 2 r (f&N^s + h.c), where t is the tunneling matrix 
element (assumed to be spatially constant). Thus, the 
effective action for the NS system consists of a sum of 
the actions of the individual systems, Sn and Ss, and 
a coupling term that in the weak tunneling limit can be 
linearized. I.e., the full action reads (see e.g. |33|) 



S 



d 2 r Str D N (dQ N 

Ds(dQs) 2 - 4(ie<r : 
^iQnQs 



&3 C Qn 



a? - Aof )Qt 



(19) 



where 7 = \t\ 2 T represents the transparency of the inter- 
face. Furthermore d — d — iA.[a^ H , .], where the appro- 
priate gauge for the vector potential A(r) will be spec- 
ified shortly. In general, the order parameter may be 
complex, A = | A| exp^xer™]. Note that A here is the 
sclf-consistently determined order parameter. The pres- 
ence of the normal layer leads to a renormalization of the 
order parameter. At weak coupling, however, the proxim- 
ity induced suppression of the order parameter is small. 
Therefore, we concentrate on the much more pronounced 
effect on the quasi-particle DoS. 

Subjecting the action ( [i"9| ) to a saddle point analysis, 
one obtains the following coupled Usadel- like equations, 



D N d(Q N dQN) 
DshQsdQs) ~ [iea p 3 



Tl il ®<T™,Q N }= 1 [Q s ,Q N ], 
f-Aaf,Qs]= 7 [Q J v,Qs]. 



As a simple guiding example, let us first consider the 
properties of the system in the absence of a magnetic 
field, where the order parameter can be chosen to be real. 
Now, in situations where the superconducting terminal is 
represented by a bulk system, the latter simply acts as a 
boundary condition for the normal region. However, in 
the present case, the single superconducting channel is 







itself heavily influenced by the contact with the normal 
region. As a result, Fominov and Feigel'man (l3) have 
shown that a gap develops in the normal region while in 
the superconductor quasi-particle states at energies down 
to the size of the proximity effect induced gap are gener- 
ated. 

To see this explicitly, let us employ the Ansatz Qx — 
cosh^xc™ <8> erg' + i sinh #x with Ox homogeneous 
(X = N,S). In this case, the saddle point equations 
reduce to 

— iesinh^jv = Tsinh^jv — 
— i(esmh8s — Acosh#s) = 7sinh(#,s — On). 

If the two systems are decoupled, 7 = 0, the solution 
for the superconductor at energies well below the gap, 
e <C A, reads s ~ iir/2. At weak coupling, setting 
9s = in/2 + $s in the low-energy regime and expand- 
ing the equations above up to linear order in $g yields 
coth^Ar = e/7 and fls = (e — i7COsh6>jv)/A. Thus, at 
small energies, the density of states in the two layers is 
given as 

f e < 7, 

^iv(e) = v n 9fi[cosh6>jv] = < e e > 7; ( 22 ) 

v s (e) = -vn %[4 S ] = ^(e). (23) 

As expected, the superconductor induces an energy gap 
in the normal region of magnitude E g ^p = 7. Fur- 
thermore, the contact with the normal region leads to 
the appearance of quasi-particle states in the supercon- 
ductor at energies down to the proximity effect induced 
gap -Egap "C A. (Note that close to the singularity at 
e = Eg&p , the approximations above are no longer valid.) 



B A 




FIG. 6. Sketch of the hybrid NS systems considered in the 
text in the planar and cylindrical geometry. 



in the first case the magnetic field effect is due only to 
the flux enclosed between the N- and S-layer, here the 
system as a whole encloses magnetic flux which leads to 
markedly different behavior. 

A. Planar geometry 

In the planar geometry, as before, the appropriate 
gauge is the London gauge A = —(Hz + c )e y . Here 
the constant cq is determined through the condition that 
the supercurrent through a cross section of the bi-layer 
vanishes fL3fl : 

[dzi(z)=j N +is=0, (24) 

where jx = nxA-x/m. To a first approximation, rijy = 
0. Thus, no supercurrent flows in the N region and, 
therefore, the supercurrent in the S region has to vanish 
as well, i.e. Ag = which implies Cq = —Hds- 

Using the same Ansatz as for the field-free case, one 
obtains 

- -^-(Hd) 2 sinh(2(9jv) - iesmhO N = 7 smh(§ N -§s), 

—i(esinh8s — A coshes) = 7 sinh(#s — On), 

where d — — ds is the distance between the two layers. 

As pointed out earlier, being a single-channel system, 
the superconductor alone does not feel the magnetic field. 
Again we are interested in the DoS at energies well below 
the gap, e -c A. As in the field free case, an expansion 
in ds = 0s — itt/2 leads to 

kjv sinh(2(?/v) + 2i(esinh0x — 7COsh(?/v) = 0, (25) 

where kn = D]y(Hd) 2 /2. Furthermore, i?s = 
(e — 17 cosh#/v)/A as before. Thus, in the two chan- 
nel case, the magnetic field leads to a suppression of the 
proximity effect. Eq. ( p5| ) shows that the closing of the 
induced gap is described by the AG theory, where the 
relevant parameter is (jv = i^n/i- Therefore, one finds 
that the characteristic field (n(H c ) = 1 causing the prox- 
imity effect induced gap to vanish is much weaker than 
the field necessary to drive the superconductor into the 
gapless phase (i.e. taking into account the finite width of 
the individual layer!). 

B. Cylindrical geometry 



How do these characteristic features of the proximity 
effect change in the presence of a magnetic field? In the 
following, we will consider two different geometries as 
depicted in Fig. ||. In Sec. VA, a plan ar NS bi-layer is 
investigated. Subsequently, in Sec. VB, we study a setup 
where the bi-layers are wrapped around a cylinder with 
the magnetic field directed along the cylinder axis. While 



While in the planar geometry a non-vanishing super- 
current is forbidden, in the cylindrical geometry a super- 
current can flow around the cylinder: in contrast to the 
previous case, the system now encloses magnetic flux. 

Starting with a single superconducting layer, for the 
cylindrical geometry, the most convenient gauge to 
choose is the symmetric gauge A = x r = —^Hre v , 
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where we use cylindrical coordinates (r, ip, z). Now the 
phase degree of freedom has to be taken into account, 
and a more general Ansatz for the matrices Q solving 
the saddle point equations is needed: 



Qs = cosh 6 S erf (g)cr!f + isinhflse" 1 ^ 3 a 



where n G Z, and the phase matches the phase of the 
order parameter, A = |A| exp[inLpa™\. Substituting this 
Ansatz into the saddle point equation, one obtains 



D s {n- 



2r§ 



sinh(20 s ) + 2i(esinh0 s -Acosh6's) = 0. 



Thus, even in the absence of the NS coupling, the mag- 
netic field affects the properties of the superconductor. 
The integer n has to be chosen such that it minimizes 



|n + Hr%\, and, therefore, by increasing the magnetic 
field, the parameter governing the gap suppression, 



(s(H;n) 



D s (n 



2r|A 



varies periodically between and £!? ax = Z?s/(8r|A) = 
(£/ds) 2 , where ds is the diameter of the cylinder. The 

size of the energy gap is E gap — A(l — Cs^ 3 ) 3 ^ an d the 
superconductor enters the gapless phase when Cs = 1. 
This condition can only be fulfilled if ds < £. 

As pointed out earlier, the magnetic field not only sup- 
presses the energy gap, but also renormalizes the order 
parameter A. I.e. the order parameter in the formulae 
above has to be determined self-consistently. At T = 0, 
the self-consistency equation can be cast in the form pi] 



ln { A^J = 1 -arcoshCs - § (Cs arcsinQ 1 - J I - Q 2 ] Cs > 1- (26) 



T /T° 

c ' c 




|A|/A 



FIG. 7. The order parameter A and the transition tem- 



perature T c as a function of the magnetic flux $ 
threading the cylinder. 



TvHrt 



Thus, the periodic modulation of £ also leads to a peri- 
odic modulation of the order parameter. Similarly, the 
transition temperature T c , which obeys 0] 



ln 



TP 



-ij> 



1 



2ttT c 



(27) 



is a periodic function of the applied field as first ob- 
served by Little and Parks Jl5[] . Furthermore, in small 
rings superconductivity is completely suppressed in a 
certain range of magnetic fields around half-integer flux 
quanta threading the cylinder pq| . Only very recently 
it has been possible to manufacture small enough cylin- 
ders, where this prediction could be verified experimen- 
tally |l7j . Fig. [7] shows the order parameter at T = and 
the transition temperature for a cylinder with ds < £0 



(where £0 is the coherence length at H — 0), plotted 
against magnetic flux $/</)o, where $ = irHr s . As ex- 
pected, both vanish around half-integer flux-quanta while 
they reach their unperturbed maximal values at integer 
flux quanta. The region around Q/(f>o — 1/2 is magnified 
in Fig. |^, where the energie gap and the order parameter 
are plotted. The system shows a cross-over S — > gapless 
S — > N — > gapless S — > S by increasing the magnetic field. 



E /A n 
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gapless 
region 



|A|/A 
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Mo 

FIG. 8. Energie gap -E gap and order paramter A in the 
vicinity of $/^o = 1/2. 



Now adding the normal layer, we again concentrate 
on energies much smaller than the gap. With the same 
approximation used earlier, the equation for the normal 
region reads 



D N (n- 



N 



2r% 



sinh(26 , Ar) + 2z(esinh#7v — 7 cosh 6^) = 0. 
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Once again, the mean-field equation assumes the form 
of an AG equation with the parameter £jv(ff;n) = 
Djsr(n + Hr%) 2 /(2rjjj). Comparing the two values Cs 
and Cn, we find that C™ ax ^ Cs^S ke. the proximity 
gap is suppressed before the superconductor itself would 
enter the gapless regime. At the same time the solution 
for the superconductor takes the form 

e — *7 cosh On 
S= A(l-Cs) ' 

yielding i>s = ^n7/A(l — Cs) _1 cosh#Ar. I.e. the com- 
bined influence of the presence of the normal region and 
the magnetic field leads to an enhanced density of states 
at low energies in the superconductor. 

This concludes our discussion of the mean-field prop- 
erties of thin disordered NS hybrid systems. Taking into 
account the influence of fluctuations, it is straightforward 
to see that the low-energy properties of the quasi-particle 
states are dictated by the same theory obtained in the 
previous section. Here, we assume the disorder in the 
N- and S-channel to be uncorrelated which violates V z - 
invariance. Thus, the gapless hybrid system is described 
by symmetry class C. 



in the DoS of the normal layer, and (ii) a finite den- 
sity of states in the superconductor at energies -Egap < 
e < A. In this geometry, a parallel magnetic field sup- 
presses the induced proximity gap £ ga p. The charac- 
teristic field H C (N) determining the occurrence of the 
gapless phase is greatly reduced as compared to the field 
H C (S) that drives the superconductor into the gapless 
phase, being roughly H C (N)/H C {S) « {E^J/A) 1 / 2 . In 
a cylindrical geometry, the energy gap shows a peri- 
odic modulation with the magnetic field reminiscent of 
the Little-Parks effect: if the cylinder encloses multi- 
ples of the flux quantum 0o> this can be compensated 
by the phase of the order parameter. Thus, the varia- 
tion of the energy gap is determined by the effective field 
H c ff — min„ e z|7J + n/r||. The gapless phase can only be 
reached in sufficiently small systems, where the diameter 
of the cylinder fulfills the relation ds < £. 
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VI. CONCLUSIONS 

To conclude, we have cast the properties of a dis- 
ordered thin superconducting film subject to a parallel 
magnetic field in the framework of a statistical field the- 
ory. In the mean-field approximation, known results from 
the AG theory || are recovered. The same phenomenol- 
ogy applies to diffusive films as well as films with colum- 
nar defects. In the diffusive case, we have shown that - 
within the gapped phase - taking into account inhomo- 
geneous instanton solutions of the saddle point equation, 
the hard gap is destroyed. By analogy with the magnetic 
impurity problem j?J , exponentially small tails within the 
gap region appear. The same is to be expected for the 
columnar defects. (For M > 2 the coupling between the 
different subbands complicates the analysis. However, 
the general behavior should not be affected qualitatively.) 

Within the gapless phase, the Berry-Robnik symme- 
try phenomenon leads to different low-energy proper- 
ties. As confirmed by numerics, for the diffusive film, 
the DoS vanishes quadratically for e — * (class C) as 
one might expect for superconducting systems where T- 
invariance is lifted. However, in the presence of only 
columnar defects, the DoS at small energies is linear in e 
(class CI) , a behavior characteristic of systems which pos- 
sess time-reversal invariance. Although the 'Pz-symmetry 
cannot prevent the gradual destruction of superconduc- 
tivity by the magnetic field, some compensation for the 
T-breaking is still effective. 

In NS bi-layers, we have shown that the coupling be- 
tween the two systems leads to (i) an energy gap E g &J 
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